%% setup
clear;clc;close all

% qq parameters
tau_q = 2;
tau_e = 2;
alpha = 0.9; 

% income parameters
y_low = 0;
y_high = 1;
p = 0.1;

%% 3 periods
T=3;
declining_fertility = 0.04;
phi = 1 - declining_fertility * (1:T-1);

simulate_A_star
f1 = fig_Astar(lambda_grid,A_star_grid,T);
print(f1,'-depsc','-painters','-noui','-r600','Astar_simulation_T3.eps')

%% 6 periods
T=6;
declining_fertility = 0.04;
phi = 1 - declining_fertility * (1:T-1);

simulate_A_star
f1 = fig_Astar(lambda_grid,A_star_grid,T);
print(f1,'-depsc','-painters','-noui','-r600','Astar_simulation_T6.eps')
